Thermo-spin effects in a quantum dot connected to ferromagnetic leads 
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We study a system composed of a quantum dot in contact with ferromagnetic leads, held at 
different temperatures. Spin analogues to the thermopower and thermoelectric figures of merit 
are defined and studied as a function of junction parameters, ft is shown that in contrast to 
bulk ferromagnets, the spin thermopower coefficient in a junction can be as large as the Seebeck 
coefficient, resulting in a large spin figure of merit. In addition, it is demonstrated that the junction 
can be tuned to supply only spin current but no charge current. We also discuss experimental 
systems where our predictions can be verified. 
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Introduction. - Thermoelectricity - the relation between 
a temperature bias and a voltage bias - is a very old 
problem of solid-state physics. It has gained renewed in- 
terest in recent years^"^ due to the prospect of utilizing 
nanostructures to develop high efficiency thermoelectric 
converters''^. Theoretical models have been put forward 
for the thermoelectric transport through quantum point 
contacts, ^'^ quantum dots,*^"'^^, molecular junctions^^ 
and other strongly correlated nanostructures^^. 

Recently, Uchida et al}^ have measured the spin equiv- 
alent of the charge Seebeck coefficient, namely the spin 
Seebeck effect, in which a temperature difference between 
the edges of a bulk ferromagnetic (FM) slab induces a 
spin- voltage difference and generates spin current. These 
authors have suggested using this effect to construct a 
spin-current source for spintronic devices^'. However, 
the spin-Seebeck coefficient in this experiment was mea- 
sured to be four orders of magnitude smaller than the 
charge Seebeck coefficient. In addition, the temperature 
difference unavoidably generates also a regular voltage 
bias across the sample, which may preclude easy appli- 
cability in spintronic devices. 

Here we study the thermo-spin effect, i.e., the spin- 
analogue to the Seebeck effect in a nanojunction com- 
posed of a quantum dot (either a molecule or a semi- 
conductor quantum dot structure) placed between two 
FM leads. The charge transport properties of such 
systems have been studied both theoretically-'^^ and 
experimentally^^' We define the spin analogues of the 
thermo-electric coefficients, and show that in this par- 
ticular case the spin- and charge-Seebeck coefficients are 
of the same order of magnitude. We also calculate the 
thermo-spin figure-of-merit (FOM), and show it to be 
relatively large, indicating high heat-to-spin- voltage con- 
version efficiency. Finally, it is demonstrated that the 
system parameters can be tuned such that a large spin- 
current can be generated without any charge current, 
thus making this system ideal for spintronic applications. 

Definitions of spin-thermal coefficients. - Consider a 
system composed of some structure (for instance a quan- 
tum dot) placed between two FM leads, which we assume 
have the same magnetization alignment. The system is 
held at a temperature Tl = Tr = T. In linear response, 



the thermo-electric Seebeck coefficient is defined as mi- 
nus the ratio between the voltage bias SV and the applied 
temperature bias ST that generates it (in the absence of 
charge current). In a spin system out of equilibrium one 
can define a spin- voltage bias as SVs = A/.tR — A/il, where 
A/ij, = /ii^i — fJ-j^i, with fXsu the electro-chemical potential 
of the spin-s on either right or left of the the quantum 
dot. We expect this bias to be essentially zero when mea- 
sured in the bulk of the FMs, but it may acquire a finite 
(albeit possibly small) value in proximity to the quantum 
dot. 

To derive the spin-Seebeck coefficient, we consider a 
system in which there is both an infinitesimal temper- 
ature bias and spin- voltage bias. The charge and spin 
currents are defined as / = /| + /|, /<; = /| — I^, re- 
spectively (note that they have the same dimensions). 
In linear response, the spin-current is given by Is — 
GsSVs -f LtST where the response coefficient Lt is re- 
lated to the fact that a temperature gradient can induce 
both a spin flow and energy flow'^'^*. Setting Is = 0, we 
find the spin-Seebeck coefficient 
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Once Ss and Gs are defined, one may define a spin-FOM, 
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where the thermal conduc- 

tance of the system, which has an electron contribution 
and a phonon contribution. The absolute value is taken 
because the spin-conductance Gs may be negative. In 
analogy with charge transport, one expects that a system 
with ZgT > 1 is a good heat-to-spin- voltage converter^. 
Model. ~ The model consists of a quantum dot between 
two FM leads. The corresponding Hamiltonian of the 
system is 

H = ^ ^ (eks - liu)cli,sC^ks + '^'^'^l^'^ 
fc,i/=L,Rs=T,i s=Ta 

+Un^ni+ J2 II (7.fc«c^fesrf^ + H.c.) ,(2) 

fe,i/=L,Rs=T,i 

where cj,^,^. creates an electron in the ^ = L, R lead with 
spin s and energy Sks (the energy depends on spin due 
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to the FM splitting) , dl creates an electron in the dot 
with spin s, fis = d\ds is the number operator, U is 
the Coulomb charging energy, Ss is the energy level in 
the dot, which is spin dependent due to a field-induced 
Zeeman splitting, AS. The latter may originate from the 
magnetic field induced by the FM leads or by an external 
field. It may also arise from the presence of spin-dipoles 
which are dynamically formed around a nanojunction^'^°. 
^uks is the coupling between the leads and the dot. This 
is the simplest system that exemplifies the physics of this 



paper but it can also be realized in experiments^^'^''. 

If the temperatures are higher than the Kondo tem- 
perature, in the sequential tunneling approximation (i.e. 
first order in ^uks) one can describe the system by using 
rate cquations^'^, which describe the populations of the 
different states in the dot. The dot can bo cither empty 
(with probability Po), populated by a spin- up electron 
(Pi), by a spin-down electron (P2) or by two electrons 
(P3). The corresponding rate equation is 
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The rates Wn^n' describe the probability per unit time 
to transfer from state n to state n' . They are evaluated 
by noting that the rate for an electron to hop onto (off) 
the dot is proportional to the probability Pvsis) to find 
an electron (hole) in the i^-lead with spin s at an energy 
e. We assume that the coupling between the leads and 
the dot is energy independent (wide band approximation) 
and for simplicity assume that the leads are symmetric 
(it is easy to show that our results, e.g., Eqs. (6) and (7) 
do not depend on junction asymmetry). We thus have 
Pusis) = o-sfvie), where /,y(e) is the Fermi distribution of 
lead V (with the corresponding temperature and chemical 
potential). The constant parameterizes both the dot- 
lead coupling and the density of states (DOS), the spin 
dependence coming from both the FM band-shift and the 
tunneling rate'^^. We assume that there is a majority of 
spin up in the ferromagnets (and that the leads magne- 
tizations are aligned), and define 7 = Q!|/a|. Thus, 7 
encodes the difference between both the DOS and the 
tunneling rates of the diffc;rent spins. 

Thus, for example, we have (setting ks = h = l,a| = 
a and = 70;) 

Wo^i = a[/((£T - Ml)/Tl) + /((£T " m)/T^R)] 



W^o^2 = a7[/((£i - Ml)/Tl) + /((e^ - Mr)/7r)] 

W^i^a = «7[/((£^ + C/-/iL)/TL) + /((£i + C/-MR)/TR)] 

W3^2 = a[(l - /((et + U- /iL)/TL)) + 

+ (l-/((eT + t/-m)/rR))] , (4) 

and similarly for the rest of the transition rates. We 
assume that phonon-induced spin-relaxation processes in 
the dot are inhibited due to the presence of FM leads, 
and are hence slower than the spin transfer time-scale 
and may be neglected. We set the chemical potentials jiL 
and /ifl as the zero of energy, and so the dot energies are 
e^^l = e T 2ij.bAB {iib is the Bohr magneton). We will 
discuss two limiting cases of small and large Zeeman field 
AB (in the sense that 2/iBAB >> a or 2/iBAP << a). 
The dot level e may be tuned, e.g., by a gate voltage. 

The steady-state solution is obtained by equating the 
right-hand side of Eq. 3 to zero. From this solution, one 
can determine the charge current, spin current and heat 
current, using the continuity equation. For the charge 
and spin currents, one has ^ = Ir — II, and ^ = 
IsR—IsL where n is the charge on the dot and m = n^—ni 
is the magnetization of the dot. Using the rate equation 
one thus obtains (setting e = h = 1) 



I. = Po (w^-^i + - Pi (Wi% - Wi%) P2 {W^% - Wt2s) - Ps {Wt2i + Wt22) 

Is. = Po {W^^i - W^^}) - A lwi% + Wi%) - P2 {W^% + W^%) - Ps {Wt2, - Wt22) , (5) 

I 



where W!^_!^^, are scattering rates of transitions between 
the dot and the u = L,K lead. Once all the currents are 
obtained, it is a matter of algebra to obtain the difl[erent 



transport coefficients using the linear response definition 
of the spin current. 

Results. - The procedure described above allows us 
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to obtain analytic expressions for all the currents and 
thermo-electric/spin coefficients. The first result is that 
in the limit of AS 0, the charge-Seebeck coefficient 5* 
is independent of 7 and is the same as was calculated in 
Rcf. fO. The spin-Seebeck coefficient Ss is found to be 
proportional to S, 



Thus, for normal leads (7 = f) we have Ss = 0, and 
for perfect FM leads (7 = or 7 = 00) the spin and 
charge coefficients are identical (up to a sign). Eq. (6) 
shows that even for a moderate value of 7 = 0.3 we have 
Ss ^ 0.55, as opposed to the bulk case where it is orders 
of magnitude smaller^'^. fn the case of large AB, the 
situation is even more interesting, since in fact Ss may 
become larger than S. In the limit of J7 ^ 00 and at 
£ = (i.e. the leads Fermi energies at the center of the 
Zeeman splitting) we find that 



S i_,exp(^) 

For a value of AB = I.5T at T = 5K, a value of 7 = 0.3 
yields ^ « 2.6. 

Let us now turn our attention to the FOM. We have 
calculated the FOM (spin and charge) numerically. For 
this we take the following parameters. The coupling be- 
tween the dot and the leads is taken as a = I0~^ meV 
(which is typical of semiconductor quantum dots). The 
charging energy is taken to be two orders of magnitude 
larger, U = 2 meV, and we take 7 = 0.2. We add to 
the thermal conductance a phonon contribution which 

is Kph = 3ko {kq = ^-^T is the quantum of ther- 
mal conductance^^'^'^), a reasonable value for nanoscale 
junctions^''. In Fig. 1 we plot the spin-FOM {ZsT, sohd 
line, blue online) and charge-FOM {ZT, dashed line, pur- 
ple online) as a function of the dot energy level e for 
two temperatures of 2K and 4K (left and right columns, 
respectively) and for AB = lO^^^T (Fig l(a-b)) and 
AB = 1.5T (Fig l(c-d)). The first value is a typical 
field produced by regular ferromagnets (e.g., iron), and 
the second corresponds to a large field splitting, which 
may be found in rare-earth ferromagnets or be induced 
by an external magnetic field. 

From Fig. 1 one can see that the behavior of the spin 
and charge FOM is similar for small Zeeman splitting, 
and that ZsT and ZT are of the same order of magni- 
tude. The situation is different for large AB, for which at 
certain energies close to e = one may obtain small ZT 
but large ZsT. This is due to the fact that the Zeeman 
splitting in that case preserves the particle- hole transport 
symmetry (the lack of which is responsible for charge 
thermo-power) but dramatically changes the transport 
properties of different spins, and hence increases the spin- 
thermopower. 

Finally, we study the system at finite currents. In the 
bulk, a temperature gradient will inevitably induce both 
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FIG. 1: (color online) Spin-FOM {ZsT, solid line, blue on- 
line) and charge-FOM {Zt, dashed line, purple online) as 
a function of e at temperatures T — 2K (left column) and 
T — 4K (right column), for two values of Zeeman splitting, 
AB = 10"^T (a-b) and AB = 1.5T (c-d) (see text for other 
numerical values). The spin- and charge-FOM are compara- 
ble in size, and at certain parameters spin efficiency may even 
exceed that of charge. 



charge and spin voltage^^, and since the spin-Seebeck 
effect is much smaller than the charge Seebeck effect, in- 
ducing large temperature biases (to generate sizeable spin 
currents) would result in even larger voltage biases. In 
the system studied here, one can instead tune the system 
parameters such that there will be a large spin current 
but vanishing charge current. 

In Fig. 2 the spin current Is (solid line) and charge cur- 
rent / (dashed line) are plotted as a function of e. Here 
the temperature T = 5K, AB ~ 1.5T and we have added 
a constant temperature gradient ST = lOK (the spin- and 
charge- voltage biases are zero). When the charge cur- 
rent vanishes (indicated by an arrow in Fig. 2) the spin 
current remains finite. The inset of Fig. 2 shows the de- 
pendence of the spin-current, evaluated by varying the 
energy e so that the charge current vanishes, as a func- 
tion of the temperature bias ST. The magnitude of the 
spin current increases with the temperature difference, 
and attains significant values for realizable temperature 
differences, until it saturates at large temperatures (note, 
however, that the saturation temperature is comparable 
to the interaction energy U, and hence one expects that 
the sequential tunneling approximation breaks down at 
these temperatures) . The finite spin current at large tem- 
perature difference stems from the fact that while the 
right lead is held at a high temperature, the temperature 
in the left lead is still low, allowing for differences in the 
tunneling rates of the different spins to be substantial. 
We stress that a situation of finite Is but vanishing / 
can not be achieved by using only a voltage bias, but a 
temperature bias is needed. 

Our results are valid even if one considers additional 
single particle levels in the dot. In the limit of infinite 
U, in fact, Eq. (6) is exact for the case of equidistant 
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FIG. 2: (color online) Spin current (solid line, blue online) 
and charge current (dashed line, purple online) as a function 
of e at a finite temperature difference &T = lOK (see text for 
other numerical parameters). The arrow indicates the energy 
at which the charge current vanishes and the spin-current 
is finite. Inset: spin-current as a function of temperature 
difference, evaluated at the energy e at which the charge- 
current vanishes. 

levels with no Zeeman splitting. In the case with Zeeman 
splitting, we have numerically estimated Ss/ S for up to 



five levels and found that even in the presence of the 
additional levels Sg/S ^ 1. 

The ability to couple a quantum dot to FM leads^^'^^ 
and to measure a local spin bias^^ have been demon- 
strated experimentally. It is thus reasonable that the 
results presented here are accessible by future experi- 
ments. Another interesting candidate for such experi- 
ments is graphene, for which both the possibility to fab- 
ricate quantum dots^^ and to bond FM leads to measure 
spin currents^^ have been demonstrated. 

We also point out that if the leads are FM, extracting 
the spin-current (or measuring the spin- voltage) has to be 
done close to the junction, at a distance shorter than the 
spin-diffusion length of the FM leads. Possible ways to 
circumvent this difficulty include the use of half-metallic 
leads (in which the spin-diffusion length should be very 
large) or to use a normal metal in contact with a thin FM 
layer for each lead, with the FM thin layers sandwiching 
the quantum dot. 
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